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Constraints on models of the Higgs boson with exotic spin and parity using the full 
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A search for particles with the same mass and couplings as those of the standard model Higgs 
boson but different spin and parity quantum numbers is presented. We test two specific non-standard 
Higgs boson hypotheses: a pseudoscalar Higgs boson with spin-parity J^ = 0“ and a graviton-like 
Higgs boson with = 2^, assuming for both a mass of 125 GeV/c^. We search for these exotic 
states produced in association with a vector boson and decaying into a bottom-antibottom quark 
pair. The vector boson is reconstructed through its decay into an electron or muon pair, or an 
electron or muon and a neutrino, or it is inferred from an imbalance in total transverse momentum. 
We use expected kinematic differences between events containing exotic Higgs bosons and those 
containing standard model Higgs bosons. The data were collected by the CDF experiment at the 
Tevatron proton-antiproton collider, operating at a center-of-mass energy of ^/s = 1.96 TeV, and 
correspond to an integrated luminosity of 9.45 fb“^. We observe no significant deviations from the 
predictions of the standard model with a Higgs boson of mass 125 GeV/c^, and set bounds on the 
possible rate of production of each exotic state. 

PACS numbers: 13.85.Rm, 14.80.Bn, 14.80.Ec 

The observation of a narrow bosonic resonance H with 
mass near 125 GeV/c^ by the ATLAS [T] and CMS [ 2 ] 

Collaborations at the Large Hadron Collider (LHC) in 
the iL ^ 77 and H ZZ decay modes, 

and the evidence of such a particle at the Tevatron, pri¬ 
marily in association with a vector boson and in decays 
to bottom-antibottom quark pairs mi, shifted the focus 
of the Higgs boson experimental program to the determi¬ 
nation of the properties of the newly discovered particle. 

The central question that needs to be addressed experi¬ 
mentally is whether only one Higgs boson is sufficient to 
explain the observed data. Specifically, the spin and par¬ 
ity of the Higgs boson should be established in order to 
determine if it plays the role predicted for it by the stan¬ 
dard model (SM) of particle physics or if it represents 
the first hint of more exotic interactions. 

The properties of the new particle observed at the LHC 
are consistent with those predicted by the SM for the 
Higgs boson. The products of cross sections and branch¬ 
ing ratios are as predicted [IllSHZl- The decays of the new 
particle to 77 , and final states, where the 

asterisk indicates an off-mass-shell Z or W vector boson, 
provide excellent samples for testing its spin and parity 
quantum numbers J and P, due to the measurable an¬ 
gular distributions of the decay products IS IS> which 
depend on the quantum numbers of the decaying parti¬ 
cle. The tests at the LHC in the bosonic decay channels 
exclude exotic states with spin and/or parity different 


from the SM prediction of = 0+ with high confidence 
level. 

At the Tevatron, the primary sensitivity to the Higgs 
boson comes from modes in which it is produced via 
its coupling to vector bosons but decays to a pair of 
fermions. While ATLAS and CMS have reported strong 
evidence for fermionic decays of the Higgs boson nniE], 
spin and parity quantum numbers have not been tested 
in these decays. As the DO Collaboration has shown [ 12 ], 
testing the spin and parity of the Higgs boson at the 
Tevatron provides independent information on the prop¬ 
erties of this particle. 

The Tevatron data can test alternative hypothe¬ 
ses in the WH, ZH production modes with H bb, by 
examining the kinematic distributions of the observable 
decay products of the vector boson and the Higgs-like 
boson [13]. Testing the spin and parity of the Higgs bo¬ 
son in H bb decays provides independent information 
on the properties of this particle. The models tested are 
described in Ref. M- For the SM case, Higgs boson as¬ 
sociated production is an S-wave process (i.e. the VH 
system is in a state with relative orbital angular momen¬ 
tum L = 0, where V = W or Z), with a cross section 
that rises proportionally to the boson speed (3 close to 
threshold. Here /3 = 2p/y/s, where p is the momentum 
of the Higgs boson in the VH reference frame and ^/s is 
the total energy of the VH system in its rest frame [14] . 
In the 0“ case, the production is a P-wave process and 
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the cross section rises proportionally to There are 
several possible = 2 + models, but for graviton-like 
models m, the production is in a D-wdive process, with 
a cross section that rises proportional to /3^. This depen¬ 
dence of the cross section on the spin-parity quantum 
numbers provides good kinematic leverage for discrimi¬ 
nating exotic from SM Higgs boson production, since the 
exotic production rate is enhanced faster than the SM 
one at larger corresponding to a larger invariant mass 
of the final state system and higher momenta of the de¬ 
cay products. The models studied predict neither the 
production cross sections for pp ^ WH^ ZH nor the de¬ 
cay branching fraction B{H ^ bh). Instead, the authors 
suggest m to purify a sample of Higgs boson candidate 
events and to study the invariant masses of the Wbb and 
Zbb systems, which differ strongly among the 0+, 0“, and 
2 + models. 

The study of the properties of a purified signal sample 
with minimal sculpting of the kinematic distributions is 
effective at the LHC in the H ZZ mode, 

which has a signal-to-background ratio s/b exceeding 2 : 1 . 
However, this is not the case for the Tevatron, where the 
SM Higgs boson searches typically have a 5/6 of 1:50 [15]. 
With the use of multivariate analyses (MVAs), small sub¬ 
sets of the data sample can be purified to achieve a 5/6 
ratio of ~ 1 : 1 . Since the events in these subsets are se¬ 
lected with MVA discriminants that are functions of the 
kinematic properties of signal and background, their dis¬ 
tributions are highly sculpted to resemble those predicted 
by the SM Higgs boson, and thus are not optimal in test¬ 
ing alternative models. 

The strategy chosen for this Letter is to generalize the 
CDF searches for the SM Higgs boson in the WH 
iubb mode [16], the ZH £~^i~bb mode [17], and the 
WH + ZH ^T^b [18] mode [19], where the Z bo¬ 
son decays into a neutrino pair or the charged lepton 
from the W-boson decay escapes detection. In the last 
case, ZH i^i~bb events may be reconstructed as ^rbb 
events if both leptons fail to meet the identification crite¬ 
ria. The generalization involves searches for pseudoscalar 
(J^ = 0 “) and graviton-like (J^ = 2 +) bosons (denoted 
X here), using MVA techniques similar to those devel¬ 
oped for the SM searches. Admixtures of SM and exotic 
Higgs particles with indistinguishable mass are also con¬ 
sidered, where exotic and SM production do not interfere 
due to different spin-parity quantum numbers. We set 
limits on the production rate times the decay branching 
ratio B{X ^ bb of the exotic boson assuming a produc¬ 
tion cross section and decay branching ratio of the exotic 
boson as predicted by the SM for the Higgs boson. We 
also test the hypotheses of the exotic models by compar¬ 
ing the data with the predictions. 

The CDF H detector is described in detail else¬ 
where EOIEI]. Silicon-strip tracking detectors [ 22 ] sur¬ 
round the interaction region and provide precise mea¬ 
surements of charged-particle trajectories in the range 


1 77 1 <2 [23]. A cylindrical drift chamber provides full 
coverage over the range \r]\ < 1. The tracking detectors 
are located within a 1.4 T superconducting solenoidal 
magnet with field oriented along the beam direction. 
The energies of individual particles and particle jets are 
measured in segmented electromagnetic and hadronic 
calorimeters arranged in a projective-tower geometry sur¬ 
rounding the solenoid. Tracking drift chambers and scin¬ 
tillation counters are located outside of the calorimeters 
to help identify muon candidates [24]. The Tevatron col¬ 
lider luminosity is measured with multicell gas Cherenkov 
detectors [25] . The data set used in the analyses reported 
in this Letter corresponds to an integrated luminosity 
of 9.45 fb“^. The data are collected using a three-level 
online event selection system (trigger). The first level, 
relying on special-purpose hardware [26] , and the second 
level, using a mixture of dedicated hardware and fast soft¬ 
ware algorithms, reduce the event accept-rate to a level 
readable by the data acquisition system. The accepted 
events are processed online at the third trigger level with 
fast reconstruction algorithms, and recorded for offline 
analysis m- 

To predict the kinematic distributions of SM Higgs bo¬ 
son events, we use the pythia [28] Monte Carlo (MC) 
program, with CTEQ5L [29] parton distribution func¬ 
tions (PDFs) of leading order (LO) in the strong cou¬ 
pling parameter ag. We scale these MC predictions to 
the highest-order cross section calculations available. To 
predict the exotic signal kinematic distributions, we use 
a modified version of madevent [30] provided by the 
authors of Ref. m- 

The predictions for the SM WH and ZH cross sec¬ 
tions m are based on the next-to-leading order (NLO) 
calculation of v2hv [32] and include next-to-next- 
to-leading order (NNLO) quantum chromodynamical 
(QCD) contributions [33], as well as one-loop electroweak 
corrections [34] . In the predictions for the decay branch¬ 
ing fractions of the SM Higgs boson [35] [36] , the partial 
decay widths for all decays except to pairs of W and Z 
bosons are computed with hdecay [37] , and the WW and 
ZZ decay widths are computed with prophecy4f [38] . 
The relevant rates are awH = (129.5 ± 9.8) fb, azH = 
(78.5 ± 5.9) fb, and B{H bb) = (57.8 ± 1.0)%. The 
uncertainties on the predicted branching ratio from un¬ 
certainties in the bottom-quark mass, and missing 
higher-order effects are estimated in Refs. [33 HQ]. 

We model SM processes and instrumental backgrounds 
using data-driven and MC methods. Simulated diboson 
( WW^ WZ^ ZZ) MC samples are normalized using the 
NLO calculations from mcem m- For tt we use a pro¬ 
duction cross section of 7.04 ±0.7 pb [42], which is based 
on a top-quark mass of 173 GeV/c^ and MSTW 2008 
NNLO PDFs [ 43 ]. The single-top-quark production cross 
section is taken to be 3.15 ± 0.31 pb [44]. The normal¬ 
ization of the %±jets and lF±jets MC samples is taken 
from ALPGEN [45] corrected for NLO effects, except in 
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the case of the WH ^ £ubb search. The normalization 
of the IT+jets MC sample in the WH ^ £vbb search, and 
the normalization of the instrumental and QCD multijet 
samples in all searches, are constrained from data sam¬ 
ples where the expected s/b ratio is several orders of mag¬ 
nitude smaller than in the search samples. The quality 
of background modeling is shown in final-state invariant 
mass distribution plots included in the Supplemental Ma¬ 
terial to this Letter, which show good agreement with the 
data in all cases. 

The analyses used to search for the exotic pseudoscalar 
and graviton-like Higgs bosons are modifications of the 
searches for the SM Higgs boson, optimized for sepa¬ 
rating the exotic signals from both the SM background 
sources and the possible SM Higgs boson signal. They 
use the most recent and efficient CDF algorithm, HO- 
BIT [46], for identifying jets from the hadronization of 
bottom quarks (6-tagging). HOBIT is a multivariate 
classifier that uses kinematic properties of reconstructed 
trajectories of charged particles (tracks) associated with 
displaced vertices, the impact parameters of the tracks, 
and other characteristics of reconstructed groups of colli¬ 
mated particles (jets) that help separate 6-jets from light- 
flavored jets. The HOBIT classifier does not perform well 
for jets with > 200 GeV and the data-based calibra¬ 
tion procedures associated with it suffer from greater un¬ 
certainties in this kinematic region. We therefore do not 
tag jets with > 200 GeV. The same tight (T) and 
loose (L) tag requirements are used as in the SM Higgs 
analyses. 

In each final state, the search channels are subdivided 
according to the number of jets, the lepton category, 
and the 6-tag category. The WX £ubb events are 
divided into 15 subchannels, corresponding to the TT, 
TL, IT, LL, and IL tagging categories of the two jets, 
for each lepton category: central leptons (electrons or 
muons), forward electrons, and isolated-track leptons. 
The ZX i~^i~bb events are divided into 16 sub¬ 
channels, corresponding to the TT, TL, IT, LL tagging 
categories in the two- and three-jet final states, sepa¬ 
rately for Z e+e“ and Z events. The 

WX^ZX ^T^b events are divided into 6 subchannels 
corresponding to the TT, TL, and IT tagging categories 
in 2-jet and 3-jet final states. A total of 37 analysis chan¬ 
nels are defined. The expected and observed event yields 
in all channels are summarized in Table [l[ summed over 
lepton, jet, and 6-tag categories. 

Two discriminant functions are defined for each sub¬ 
channel, one to separate the exotic Higgs boson sig¬ 
nal (separately defined for the 0“ and the 2+ signals) 
from the backgrounds, and the other as the discrimi¬ 
nant used in the search for the SM Higgs boson. For 
the ZX i~^£~bb analysis, only the exotic discriminant 
is used. The exotic signal discriminants have either 
(the invariant mass of the final-state system) among their 
input variables or Ht (the sum of all transverse energies 


reconstructed in the final state, including muon energies 
and ^t)- Distributions of the discriminant functions for 
all search channels are shown for the data and simula¬ 
tion in the Supplemental Material to this Letter. Since 
the events are primarily classified to test for the exotic 
models, the SM Higgs interpretation of the data will not 
be the same as in the searches optimized for the SM Higgs 
boson. 

TABLE I: Expected and observed event yields for all channels. 
The difference between the 0“ and 2^ exotic yields is due to 
different signal acceptances. 


Process 

£+£-bb 

£vbb 

^j^bb 

V + Xo- 

8±1 

49±4 

81±6 

V + X 2 + 

7±1 

43±4 

65±5 

VH 

7±1 

33±3 

40±3 

U-l-jets 

820T141 

23323T2860 

9193T2273 

Dibosons 

72T11 

1288T148 

544T66 

Top 

222 T22 

2053T211 

1935T164 

QCD 

58T21 

2406T603 

16283T1447 

Total bkg 1172±272 29070±3037 27956±3188 

Observed 

1182 

26337 

28518 


To summarize the data in the large number of con¬ 
tributing channels, we follow Ref. [4]. We sum the con¬ 
tents of bins with similar s/b ratios over the output his¬ 
tograms of all channels. Fig. shows the comparison of 
the data with the best-fit background predictions and the 
summed signals, separately for the SM Higgs and exotic 
boson signals. The signal strength modifier is denoted by 
Mexotic, which multiplies the SM signal strength to pre¬ 
dict the rate in the exotic model under test. Both distri¬ 
butions show agreement between the background predic¬ 
tions and the observed data over five orders of magnitude. 
No evidence for an excess of exotic signal-like candidates 
is seen. 

A number of systematic uncertainties among the vari¬ 
ous analyses affect the sensitivity of the final result. All 
correlations within and between channels are taken into 
account in deriving the following combined limits, cross 
sections, and p-values. Uncertainties of 5% jSU |47| on 
the inclusive WH and ZH production rates are estimated 
by varying the factorization and renormalization scales. 
We assign uncertainties to the Higgs boson branching ra¬ 
tios as calculated in Ref. m- Since the exotic signals 
are normalized to the SM Higgs cross section, the same 
relative uncertainties are assumed for the exotic produc¬ 
tion. The largest sources of uncertainty on the dominant 
backgrounds are the rates of Ufheavy-flavor jets. The 
resulting uncertainties are up to 8% of the predicted val¬ 
ues. Because the various analyses use different methods 
to obtain the U+heavy-flavor predictions, we treat their 
uncertainties as uncorrelated between the £ubb^ the ^t66 , 
and £~^£~bb channels. We use simulated events to study 
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FIG. 1: Distributions of logiQ(s/6) for the data from all 
contributing Higgs boson search channels for a boson mass 
of 125 GeV/c^ for (a) the 0“ search and (b) the 2^ search. 
The observed numbers of events are represented by the points, 
and the expected exotic signals are shown as histograms with 
Mexotic = 1 stacked on top of the backgrounds, which are ht to 
the data within their systematic uncertainties. The expected 
s/b ratios of the exotic signal over the background yield are 
used to rank analysis bins. The background predictions do not 
include the contributions from the SM Higgs boson, which 
are shown as separate histograms, not stacked. The error 
bars shown on the data correspond to the square root of the 
observed data count. Underflows and overflows are collected 
into the leftmost and rightmost bins, respectively. 


the impact of the jet-energy-scale uncertainty [48] on the 
rates and shapes of the signal and background expecta¬ 
tions. We treat the jet-energy-scale variations uncorre¬ 
lated among the three analyses in the combined search. 
Uncertainties on lepton identification and trigger efficien¬ 
cies range from 2% to 6% and are applied to both signal 
and MC-based background predictions. The uncertainty 
on the integrated luminosity is 6%, of which 4.4% orig¬ 
inates from detector acceptance uncertainties and 4.0% 
is due to the uncertainty on the inelastic pp cross sec¬ 
tion m- The luminosity uncertainty is correlated be¬ 
tween the signal and MC-based background predictions. 

Bayesian exclusion limits at 95% credibility level 


(c.L.) m on the production rates times the branch¬ 
ing fraction B{X bh) for 0“ and 2+ Higgs bosons are 
reported in Table [H] both separately for each channel 
and combined, in units of the SM Higgs boson produc¬ 
tion rate. The limits are computed from a likelihood 
defined as the product of the probability densities for 
the bin contents of the MVA histograms over all bins 
of each histogram and all channel histograms, assuming 
Poisson probability densities for the bin contents, uni¬ 
form prior densities for the SM and exotic signal strength 
modifiers /iexotic and /isM, and Gaussian prior densities 
for the nuisance parameters describing systematic uncer¬ 
tainties. Posterior densities and upper limits on the SM 
and exotic Higgs boson rates are obtained from pseudo¬ 
experiments (PEs), where in each PE the likelihood is 
integrated over the nuisance parameters and then it is 
maximized. The medians of the distributions of results 
from PEs are used as the most probable values. The SM 
ratio between WH and ZH production rates is assumed 
when combining WX and ZX searches. Limits are listed 
either assuming that the SM Higgs boson is present as a 
background, or absent. Since the exotic 0“ and 2+ sig¬ 
nals populate kinematic regions different from those of 
the SM Higgs boson, and since the SM Higgs boson pro¬ 
duction rate is small, the expected and observed limits 
on the exotic rates are very similar whether the SM Higgs 
boson is present or not. The observed combined limits 
are somewhat stronger than expected, with an exclusion 
rate of /iexotic < 0.32 in the 0“ case (approximately a 
one standard deviation deficit), and /iexotic < 0.35 in the 
2 + case (approximately a two standard deviation deficit). 
The channel carries the largest weight in the com¬ 

bination. A number of candidates somewhat lower than 
expected appear in the most signal-like bins of the exotic 
discriminants in this channel. The two-dimensional cross 
section fits, which allow for arbitrary rates of both SM 
and exotic Higgs bosons to be simultaneously present, are 
shown in Eig. separately for the 0“ and 2+ searches. 

TABLE H: Limits at 95% G.L. on 0“ and 2^ boson production 
assuming no SM Higgs boson background. In parentheses are 
the limits assuming SM Higgs boson background. 



0 


2 

+ 

Channel 

Obs 

Median exp 

Obs 

Median exp 


[limit/ Hsm] 

[limit/ Hsm] 

[limit/ Hsm] 

[limit/ Hsm] 

iubb 

0.59 (0.55) 

0.74 (0.78) 

1.05 (0.99) 

1.01 (1.03) 

i^i-bb 

1.86 (1.77) 

1.46 (1.52) 

1.57 (1.49) 

1.59 (1.61) 

^Tbb 

0.49 (0.43) 

0.68 (0.69) 

0.41 (0.37) 

0.79 (0.83) 

Combined 

0.32 (0.28) 

0.44 (0.45) 

0.35 (0.31) 

0.54 (0.56) 


We report the observed values and the expected dis¬ 
tributions of the log-likelihood ratio (LLR) [50] in the 
SM and the exotic hypotheses and list the combined 
results in Table m The Table includes the p-values 
for the null and test hypotheses, defined as the condi¬ 
tional probabilities Pnuii = P(LLR < LLRobs|SM) and 
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production, and assuming only exotic production. 

Acknowledgments 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu¬ 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
0 1 2 3 4 5 6 Council of Canada; the National Science Council of the 

Republic of China; the Swiss National Science Founda¬ 
tion; the A.P. Sloan Foundation; the Bundesministerium 
fiir Bildung und Forschung, Germany; the Korean World 
Class University Program, the National Research Foun¬ 
dation of Korea; the Science and Technology Facilities 
Council and the Royal Society, United Kingdom; the 
Russian Foundation for Basic Research; the Ministerio de 
Ciencia e Innovacion, and Programa Consolider-Ingenio 
2010, Spain; the Slovak R&D Agency; the Academy of 
Finland; the Australian Research Council (ARC); and 
the EU community Marie Curie Fellowship Contract No. 
302103. 

0 1 2 3 4 5 6 

'^SM - 


0.5 

0.4 

0.3 

0.2 

0.1 


' ' ' ' I ' ' ' J M 

• Best Fit 
aSM 

■ 68% CL 
□ 95% CL 




(b) 


0.5 r 
0.4 r 
0.3 r 
0.2 
0.1 


m Best Fit 
aSM 

068% CL 
□ 95% CL 


(a) 


FIG. 2: Combined two-dimensional posterior density of the 
measured (a) 0“-vs.-0^ and (b) 2^-vs.-0^ cross sections, nor¬ 
malized to the SM predictions. 


Ptest = P(LLR > LLRobs I exotic), respectively, the val¬ 
ues of CLs = Ptest/(1 — Pnuii), and the equivalent number 
of Gaussian standard deviations 2 ; corresponding to each 
p-value, defined by p = [1 — erf(z/v^)]/2 [50]. There is a 
deficit in the observed number of events in the signal-like 
bins of the exotic discriminant, which is visible in Fig. 
in both the 0“ and the 2+ searches. The dominant contri¬ 
bution to this deficit comes from the WX ZX 
search. This deficit in the exotic search is not evidence 
against the SM Higgs boson, as the exotic search tests for 
events with different kinematic properties (high 
than those of the SM Higgs boson. Indeed, the combined 
cross section fit, shown in Fig. is consistent with the 
SM Higgs boson rate with a discrepancy of less than 0.5 
standard deviations. 

In conclusion, we search in the entire GDF data sam¬ 
ple for Higgs-boson-like particles of the same mass, pro¬ 
duction and decay modes, and production rates as the 
discovered SM Higgs boson candidate, but with 0“ or 2+ 
spin-parity quantum numbers. We observe no significant 
deviations from the SM predictions with a Higgs boson 
of mass tuh ~ 125 GeV/c^, and set bounds on the pos¬ 
sible rate of production of 0“ and 2+ exotic states, both 
allowing for an admixture of SM production and exotic 
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TABLE III: LLR and p-values for the test hypotheses. The SM hypothesis includes a SM Higgs boson. The significances 
corresponding to the p-values and CLg are given in parentheses in units of standard deviation (s.d.). The negative signal 
significance Pnuii reflects the dehcit of signal-like events compared with the background prediction. 
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FIG. 3: ^rhb two-dimensional posterior density of the measured (a) 0 -vs.-O^ and (b) 2^-vs.-0^ cross sections, normalized to 
the SM predictions. 




FIG. 4: liybb two-dimensional posterior density of the measured (a) 0 -vs.-O^ and (b) 2^-vs.-0^ cross sections, normalized to 
the SM predictions. 
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FIG. 5: llhb two-dimensional posterior density of the measured (a) 0 -vs.-O^ and (b) 2^-vs.-0^ cross sections, normalized to 
the SM predictions. 



FIG. 6: ^rbh LLR distributions for (a) 0 and (b) 2+ hypotheses, assuming /x = 1. 



FIG. 7: lybb LLR distributions for (a) 0 and (b) 2^ hypotheses, assuming /x = 1. 
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FIG. 8: llhh LLR distributions for (a) 0 


and (b) 2^ hypotheses, assuming ji — 1. 



FIG. 9: Gombined LLR distributions for (a) 0 and (b) 2^ hypotheses, assuming /i = 1. 
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FIG. 10: Background-subtracted distribution of the discriminant histograms, summed for bins with similar signal-to-background 
ratio {s/h) over all contributing Higgs boson search channels, for mn = 125 GeV/c^, for (a) the 0“ search and (b) the 2+ 
search. The background is fit to the data in each case, and the uncertainty on the background, shown with dashed lines, is 
after the fit. The exotic signal model, scaled to the SM Higgs boson expectation, is shown with a filled histogram. The SM 
Higgs boson expectation is also shown with a solid line. The error bars shown on the data points correspond in each bin to the 
square root of the sum of the expected signal and background yields. Underflows and overflows are collected into the leftmost 
and rightmost bins, respectively. 
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FIG. 11: Exclusion limits on the fraction of exotic boson in the total signal for the exotic + SM Higgs boson admixture 
hypotheses. The limits are shown as functions of the total signal production rate relative to the SM prediction for the Higgs 
boson rate. 
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FIG. 12: Discriminant function for (a) the tight-tight and (b) the tight-loose ^-jet categories of the iiybb search channel for a 
0“ state. Events are first categorized according to their exotic discriminant outputs. If these pass below a threshold of 0.5, 
then the SM discriminant function is used instead, scaled to the range [0,0.5]. In this way, kinematic regions that are the most 
sensitive to testing for the presence of an exotic Higgs boson are used first, and non-exotic-like events are used to compare the 
data with the SM Higgs boson hypothesis. 



FIG. 13: Discriminant function for (a) the loose-loose and (b) the 1-tight (right) &-jet categories of the iubb search channel for 
a 0“ state. 



FIG. 14: Discriminant function for the 1-loose 5-jet category of the iubb search channel for a 2^ state. 
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FIG. 15: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 2-jet bin of the Uhh search 
channel for a 0“ state. 
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FIG. 16: Discriminant function for (a) the loose-loose and (b) the 1-tight 5-jet categories in the 2-jet bin of the llhh search 
channel for a 0~ state. 
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FIG. 17: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 3-jet bin of the llhh search 
channel for a 0~ state. 
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FIG. 18: Discriminant function for (a) the loose-loose and (b) the 1-tight ^-jet categories in the 3-jet bin of the Mbb search 
channel for a 0“ state. 
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FIG. 19: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 2-jet bin of the prbb 
search channel for a 0“ state. Events are first categorized according to their exotic discriminant outputs. If these pass below a 
threshold (typically 0.5 on a scale of 0 to 1), then the SM discriminant function is used instead, scaled to the range [0,0.5]. In 
this way, kinematic regions that are the most sensitive to testing for the presence of an exotic Higgs boson are used first, and 
non-exotic-like events are used to compare the data with the SM Higgs boson hypothesis. 
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FIG. 20: Discriminant function for the 1-tight 5-jet category in the 2-jet bin of the ^rbb search channel for a 0 state. 
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FIG. 21: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 3-jet bin of the ^rbb search 
channel for a 0“ state. 



FIG. 22: Discriminant function for the 1-tight 5-jet category in the 3-jet bin of the ^rbb search channel for a 0 state. 
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FIG. 23: Discriminant function for (a) the tight-tight and (b) the tight-loose ^-jet categories of the iiybb search channel for a 
2^ state. Events are first categorized according to their exotic discriminant outputs. If these pass below a threshold of 0.5, 
then the SM discriminant function is used instead, scaled to the range [0,0.5]. In this way, kinematic regions that are the most 
sensitive to testing for the presence of an exotic Higgs boson are used first, and non-exotic-like events are used to compare the 
data with the SM Higgs boson hypothesis. 
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FIG. 24: Discriminant function for (a) the loose-loose and (b) the 1-tight 5-jet categories of the iubb search channel for a 2^ 
state. 
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FIG. 25: Discriminant function for the 1-loose 5-jet category of the iubb search channel for a 2^ state. 
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FIG. 26: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 2-jet bin of the Uhh search 
channel for a 2^ state. 
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FIG. 27: Discriminant function for (a) the loose-loose and (b) the 1-tight 5-jet categories in the 2-jet bin of the llhh search 
channel for a 2^ state. 
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FIG. 28: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 3-jet bin of the llhh search 
channel for a 2^ state. 
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FIG. 29: Discriminant function for (a) the loose-loose and (b) the 1-tight (right) 5-jet categories in the 3-jet bin of the iibb 
search channel for a 2^ state. 
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FIG. 30: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 2-jet bin of the prbb 
search channel for a 2^ state. Events are first categorized according to their exotic discriminant outputs. If these pass below a 
threshold of 0.6, then the SM discriminant function is used instead, scaled to the range [0,0.5]. In this way, kinematic regions 
that are the most sensitive to testing for the presence of an exotic Higgs boson are used first, and non-exotic-like events are 
used to compare the data with the SM Higgs boson hypothesis. 
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FIG. 31: Discriminant function for the 1-tight 5-jet category in the 2-jet bin of the ^t 55 search channel for a 2^ state. 
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FIG. 32: Discriminant function for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 3-jet bin of the ^rbb search 
channel for a 2^ state. 
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FIG. 33: Discriminant function for the 1-tight 5-jet category in the 3-jet bin of the ^rbb search channel for a 2^ state. 
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FIG. 34: Reconstructed Vbb transverse invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories 
of the iiybb search channel. 




FIG. 35: Reconstructed Vbb transverse invariant mass distribution for (a) the loose-loose and (b) the 1-tight 5-jet categories of 
the iiybb search channel. 



FIG. 36: Reconstructed Vbb transverse invariant mass distribution for the 1-loose 5-jet category of the iubb search channel. 













































26 




FIG. 37: Reconstructed Vbb invariant mass distribution for (a) the tight-tight and (b) the tight-loose (right) 6-jet categories in 
the 2-jet bin of the Mbb search channel. 




FIG. 38: Reconstructed Vbb invariant mass distribution for (a) the loose-loose and (b) the 1-tight 6-jet categories in the 2-jet 
bin of the iibb search channel. 




FIG. 39: Reconstructed Vbb invariant mass distribution for (a) the tight-tight (left) and (b) the tight-loose 6-jet categories in 
the 3-jet bin of the Mbb search channel. 
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FIG. 40: Reconstructed Vbb invariant mass distribution for (a) the loose-loose (left) and (b) the 1-tight 5-jet categories in the 
3-jet bin of the Mbb search channel. 




FIG. 41: Reconstructed Vbb transverse invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories 
in the 2-jet bin of the ^rbb search channel. 



FIG. 42: Reconstructed Vbb transverse invariant mass distribution for the 1-tight 5-jet category in the 2-jet bin of the ^t 55 
search channel. 
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FIG. 43: Reconstructed Vbb transverse invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories 
in the 3-jet bin of the ^rbb search channel. 



FIG. 44: Reconstructed Vbb transverse invariant mass distribution for the 1-tight 5-jet category in the 3-jet bin of the ^rbb 
search channel. 
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FIG. 45: Reconstructed bb invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories of the iubb 
search channel. Overflows are added to the content of the uppermost bin. 
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FIG. 46: Reconstructed bb invariant mass distribution for (a) the loose-loose and (b) the 1-tight 5-jet categories of the iubb 
search channel. Overflows are added to the content of the uppermost bin. 



FIG. 47: Reconstructed bb invariant mass distribution for the 1-loose 5-jet category of the iubb search channel. Overflows are 
added to the content of the uppermost bin. 
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FIG. 48: Reconstructed bb invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 2-jet 
bin of the iibb search channel. 



FIG. 49: Reconstructed bb invariant mass distribution for (a) the loose-loose and (b) the 1-tight 5-jet categories in the 2-jet bin 
of the iibb search channel. 



FIG. 50: Reconstructed bb invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 3-jet 
bin of the iibb search channel. 
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FIG. 51: Reconstructed bb invariant mass distribution for (a) the loose-loose and (b) the 1-tight 5-jet categories in the 3-jet bin 
of the iibb search channel. 



FIG. 52: Reconstructed bb invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 2-jet 
bin of the ^rbb search channel. 
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FIG. 53: Reconstructed bb invariant mass distribution for the 1-tight 5-jet category in the 2-jet bin of the ^rbb search channel. 
























































































32 



FIG. 54: Reconstructed bb invariant mass distribution for (a) the tight-tight and (b) the tight-loose 5-jet categories in the 3-jet 
bin of the ^rbi search channel. 



FIG. 55: Reconstructed bb invariant mass distribution for the 1-tight 5-jet category in the 3-jet bin of the ^rbb search channel. 











































